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Measuring Cluster Temperature Profiles with XMM/EPIC 
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Abstract. Using the PV observation of A1795, we illustrate the capability of XMM-EPIC to measure cluster 
temperature profiles, a key ingredient for the determination of cluster mass profiles through the equation of hy- 
drostatic equilibrium. We develop a methodology for spatially resolved spectroscopy of extended sources, adapted 
to XMM background and vignetting characteristics. The effect of the particle induced background is discussed. A 
simple unbiased method is proposed to correct for vignetting effects, in which every photon is weighted according 
to its energy and location on the detector. We were able to derive the temperature profile of A1795 up to 0.4 
times the virial radius. A significant and spatially resolved drop in temperature towards the center (r < 200 kpc) 
is observed, which corresponds to the cooling flow region of the cluster. Beyond that region, the temperature is 
constant with no indication of a fall-off at large radii out to 1.2 Mpc. 
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The determination of the temperature profiles of the hot 
gas in galaxy clusters is essential to derive the gas en- 
tropy distribution and to measure the total mass content, 
through the equation of hydrostatic equilibrium. These 
are key quantities for cosmological studies based on the 
properties of galaxy clusters. The baryonic mass fraction 
in clusters is a strong constraint on the density parameter 
of the Universe (Briel, Henry & Bohringer 1992; White 



This is an independent and powerful probe of f2, as well 
as, the index and normalization of the power spectrum 
of primeval fluctuations (e.g. Oukbir & Blanchard 1997). 



et al. 1993). If the mass-temperature relation is well cal- 



ibrated, the mass distribution function and its evolution 
can be derived from the observed temperature function. 



to: 



M. 



Arnaud, 



Send offprint requests 
naud@discovery.saclay.cea.fr 

* Based on observations obtained with XMM-Newton, an 
ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and the USA (NASA) 

** EPIC was developed by the EPIC Consortium led by 
the Principal Investigator, Dr. M. J. L. Turner. The con- 
sortium comprises the following Institutes: University of 
Leicester, University of Birmingham, (UK); CEA/Saclay, IAS 
Orsay, CESR Toulouse, (France); IAAP Tuebingen, MPE 
G arching, (Germany); IFC Milan, ITESRE Bologna, IAUP 
Palermo, Italy. EPIC is funded by: PPARC, CEA, CNES, DLR 
and ASI 



Furthermore, the physics of structure formation and evo- 
lution can be assessed. For example, the total mass pro- 
file gives information on the physics of gravitational col- 
lapse. Current numerical simulations under a cold dark 
matter scenario predict a universal dark matter profile 
(Navarro, Frenk & White 1997 ). The gas entropy distri- 
bution, the relative distribution of the gas and dark matter 
and departures from self-similarity constrain the thermo- 
dynamical history of the hot gas and the importance of 
additional non-gravitational effects, such as galaxy feed- 
back (e.g. Ponman, Cannon & Navarro 1999] ). 

While gas density profiles have been determined with 
good accuracy by ROSAT, the shape of the temperature 
profiles in clusters, as measured by ASCA and SAX, is 
still a matter of debate (Markevitch et al. |1998| ; Irwin et 



al. 



19 99|; W hite pOOCf Irwin & Brcgman pOOO 
20001 



Molendi 



et al. pOOq Ettori et al. gOOOJ ). It is at present unclear 
if cluster atmospheres are essentially isothermal or if the 
temperature decreases with radius (as expected from nu- 
merical simulations of cluster formation). As a result, total 
mass estimates are ty picall y uncertain by a factor of two 
(Neumann & Arnaud 1999; Horner, Mushotzky & Scharf 
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1999) and the mass profiles are too uncertain to provide 
any real constraint on theoretical models. 

Significant progress was expected with XMM- 



Newton/EPIC (Jansen et al. |200l[ Turner et al. |2001[ ), 
which has a much better sensitivity than ASCA and 
SAX and does not suffer from the large energy depen- 
dent PSF of ASCA, a major source of systematic un- 
certainty. In this paper, we illustrate the capability of 
XMM-Newton to measure cluster temperature profiles, us- 
ing the PV observation of A1795. A1795 is a bright clus- 
ter (S x [2 - 10] keV > 5 x 10- 11 ergs/s/cm 2 . Its redshift 
(z = 0.063) allows a significant coverage of the cluster 
by the XMM-Newton field of view (15' in radius which 
corresponds to 1.48 h^ Mpc or about 0.5 times the clus- 
ter virial radius). The paper is organized as follows: af- 
ter this introduction, Section 2 presents a methodology 
for spectro-imagery adapted to XMM background and 
vignetting characteristics. Our results in term of cluster 
global properties and temperature profiles are presented in 
Section 3, with particular emphasis on the various sources 
of errors. Section 4 contains our conclusions. 



2. Data Analysis 

A1795 was observed for ~ 50 ksec in Revolution 100 in 
Full Frame Mode with the EPIC/MOS and pn camera. 
We focus here on the MOS data. We generated calibrated 
event files with SASv4.1, except for the gain correction. 
Correct PI channels are obtained by interpolating gain 
values obtained from the observations of the on-board cal- 
ibration source closest in time to A1795 observation. Data 
were also checked to remove any remaining bright pixels. 
Spectra in various regions were extracted to study tem- 
perature variations. Only events corresponding to detector 
regions in view of the sky are considered (using SASv4.1). 
The region corresponding to the bright point source in the 
south was also excluded. A proper treatment of the back- 
ground and vignetting effects are essential when analyzing 
extended sources like cluster of galaxies. This is discussed 
in the following two sections. 
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Fig. 2. Combined EPIC/MOS1&2 image of A 1795 in the 
[0.3 — 10] keV energy band (logarithmic intensity). 

2.1. Background 

The XMM background shows considerable variations with 
time, with flares of various durations and intensities. As an 
example, we show in Fig. |l| the count rate in the [10 — 12] 
keV energy band, where the emission is dominated by 
the background, due to the small effective area of XMM- 
Newton/EPIC/MOS in this band. It is essential to remove 
these periods of high background, which are induced by 
soft protons from solar flares collected by the telescopes. 
Not only is the S/N highly degraded, especially at high 
energy (where the data are crucial for kT measurements), 
but it is impossible to properly subtract the background, 
since its spectrum is varying with time. We thus removed 
all frames corresponding to a count rate greater than 15 
ct/lOOs in the [10 — 12] keV band. After this data clean- 
ing, the useful observing time for MOS1 and MOS2 is 
32.1 ksec and 32.3 ksec, respectively. The corresponding 
EPIC/MOS1&2 image in the [0.3-10] keV energy band is 
shown on Fig.|^. The remaining background is dominated 
by the cosmic X-ray background at low energy and the 
Cosmic Ray (CR) induced background at energies above 
typically 1.5 keV. 

To estimate the background we used the Lockman Hole 
(LH) observations made with the thin filter (Rev 70 and 
Rev 73 for MOS1 and Rev 71 for MOS2). Data were fil- 
tered for high background and regions corresponding to 
the 10 brightest sources were excluded^]. 

It is known that the CR induced background changes 
slightly in the FOV. It is thus better to consider the same 
extraction regions in detector coordinates for the source 
and the background. Figure |^ shows the raw cluster spec- 
trum in the 4.5' — 5' region, the background spectrum 
derived from the LH observation and the corresponding 

1 Background event files, combining several high galactic lat- 
itude pointings, and thus allowing a statistically better esti- 
mate of the CR induced background, were generated after the 
completion of this work (D.Lumb, 2000) and was used in our 
analysis of Coma (this issue). The derived background is fully 
consistent with the background used in this analysis 
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Fig. 3. Summed EPIC/MOS1 and MOS2 spectra in the 
4.5' — 5' region. Green points: raw spectrum. Red points: the 
background spectrum derived from the LH observation. Black 
points: the corresponding background subtracted cluster spec- 
trum. 



background subtracted spectrum. Note the Al and Si flu- 
orescence lines in the background spectrum and the very 
hard continuum above 2 keV typical of CR induced back- 
ground. As a result, even for a bright cluster like A1795 
and relatively close to the center, the background domi- 
nates the emission at high energies (here above 7 keV). 
This is an important limitation for temperature estima- 
tion. 

The CR induced background changes slightly with 
time: the average background level obtained after flare 
cleaning varies typically by less than 10% from observation 
to observation. We considered two criteria for background 
normalization (with each camera treated independently): 
the count rate in the [10—12] keV band in the FOV (where 
sky X-ray emission is negligible) and the count rate in the 
[0.3 — 12] keV band in the region of the detector masked 
by the proton shield (no sky emission). The derived nor- 
malization factor between the LH and A1795 observations 
varies from 0.94 to 1.03 depending on the camera and 
criteria considered. In the following we use a nominal a 
normalization factor of unity. It will be varied by ±5% 
to assess systematic uncertainties due to background sub- 
traction. 

Note that this treatment allows a proper estimate of 
the CR induced background, but not of the soft X-ray dif- 
fuse background (which varies with position on the sky) 
and of the extra-galactic background (which depends on 
the absorbing hydrogen column density along the line of 
sight). However, the source usually dominates the back- 
ground except at high energy (Fig ||) and we are mostly 
sensitive to CR induced background. Furthermore, Nh 
values are similar in the directions toward A1795 and LH 
an d we check ed from ROSAT survey maps (Snowden et 
al. 1995 , 1997 ) that the X-ray emission around the LH and 



A1795 in the soft band ([0.4-1.2] keV) differ by less than 
30% on average. 

2.2. Correction for vignetting effects 
2.2.1. Methodology 

While the XMM-Newton/MOS spectral resolution does 
not show spatial variations, the effective area at a given 
energy does depend on position. This effect has to be taken 
into account in modeling spectra extracted from various 
regions. The exact method is to perform a convolution of 
a source model with the instrument response (i.e. point 
by point) and to compare the resulting modeled spectrum 
with the observed spectrum of the chosen region. This 
method is complex to implement, especially if the region 
under consideration is not assumed to be isothermal. In 
any case, it requires an a priori model for the source spatial 
distribution. 

A simpler and widely used method is to use the emis- 
sion weighted effective area in the region, derived from 
the observed global photon spatial distribution (the same 
weighting is done at each energy). The incident source 
spectrum is simply compared with the observed spectrum 
using this effective area (ARF). This method has severe 
drawbacks. First, it introduces extra noise, since the ob- 
served photon distribution used in the ARF average is 
noisy. Taking this error into account is not trivial. Even 
more important, the method introduces a bias in the de- 
termination of the spectral parameters. Since the effective 
area decreases with distance from the center, the detected 
distribution of photons (at a given energy or integrated 
over the energy band) is always more pronounced towards 
the center than the source distribution (the well known 
distortion of the observed image with respect to the in- 
cident image). If the average of the effective area is done 
using the observed distribution, too much weight is given 
to the central regions. The overall effective area is thus 
overestimated. Moreover since the variation with energy 
of the effective area depends on position (the effective area 
decreases more rapidly with distance for high energy pho- 
tons than for low energy ones), a bias is introduced. The 
overestimate of the effective area is higher at high energy 
than at low energy. The resulting effect is that the fit for 
the temperature is biased to lower values. 

We propose instead the following method, which is 
easy to implement, does not introduce bias and does not 
require any a priori assumptions about the spatial varia- 
tion of the source (either in intensity or in spectroscopic 
properties). When extracting the spectrum of a region 
Reg, we weight each photon with energy Ej falling at 
position (xj , y.j ) by the ratio of the effective area at that 
position A Xj:Vj (Ei) to the central effective area Aq^Ei). 
We can thus define a 'corrected' spectrum C(I), where C 
is the corrected photon counts in channel I centered on E\ 
with width AEj: 



C{I) = ^wj; (xj,yj) e Reg; 



(1) 
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Ei - AEi/2 < Ej < Ei + AEt/2 



Aq, (Ei) 
{Ex) 



A 



(2) 



This weighting factor has to be taken into account in 
the error estimate. The variance on C(J) is 



W i 



(3) 



Note that this method corresponds to CORRECT 
mode in EXSAS for ROSAT data reduction (Zimmermann 
et al. 1998) ) . This 'corrected' spectrum is an estimate of 
the spectrum one would get if the detector were flat. It can 
thus be compared with the total spectrum of the source 
given by the model, convolved with the instrument re- 
sponse at the center of the detector. The method is exact 
for an instrument with perfect spatial and spectral (AE) 
resolution, i.e., the energy and position of the detected 
photons correspond exactly to those of the emitted pho- 
tons. However, the method cannot introduce significant 
bias as long as the vignetting effect is small at the scale 
of the PSF and remains the same for energies differing by 
about AE. This is certainly the case for XMM-Newton. Its 
only drawback, as compared to the direct exact method, is 
a degradation of the statistical quality: the relative error 

is increased by a factor ^JJwf) / (wj) where the brackets 

denote the average over the photons in a given energy bin 
I. This is not a problem if the effective area does not vary 

much in the region being considered {\J {wj) / ' {wj) ~ 1)- 

For consistency, the background spectra were obtained 
using the same correction method as for the source. The 
CR- induced background component is not vignetted, but, 
since we extract the background and source spectra from 
the same region in detector coordinates, the correction 
factor is the same and thus cancels. 



2.2.2. Vignetting data 

For the vignetting function of the telescope, we used the 
data available on the XMM-Newton home page in August 
2000, based on a simplified model of the mirrors. Data 
in the SASv4.1 CCF are wrong and correct values de- 
rived from the last model implemented in SciSim were 
only made available after the completion of this work. We 
add the vignetting effect due to the obscuration by the 
Reflection Grating Array (RGA), as derived from SciSim 
(P.Gondoin private communication). This effect can be 
supposed to be energy independent but introduces an az- 
imuthal dependence in the vignetting function. The over- 
all vignetting is in reasonable agreement with values de- 
rived from a comparison of on- and off-axis (10' from 
centre) observations of the supernova remnant G21.5— 0.9 
(Neumann, 2000). 




channel energy (keV) 

Fig. 4. XMM and ROSAT spectra of the cluster from 
the region covering radii of 2' — 10'. Green (blue) points: 
EPIC/MOSl(2) data. Black points: PSPC data. The EPIC 
spectra are background subtracted and corrected for vignetting 
as described in Sec. ^. The red line is the best fit isothermal 
model with kT = 5.95 keV, an abundance of 0.27 times the 
solar value and a relative normalization between the EPIC and 
PSPC data of 1.05. 

2.3. Spectral fitting 

The source spectra (with errors) computed as described 
above, are binned so that the S/N ratio is greater than 3 
er in each bin after background subtraction. The spectra 
are fi tted with XSPEC using a MEKAL model (Me we 
et al. |1985| , [1986| Kaastra $992} Liedahl et al. |1995| ). Since 
the spectra are 'corrected' we can use the on axis response 
file. Note that the response matrix we used (version v3.15) 
assumes no energy dependence for the RGA transmission. 
Only data above 0.3 keV are considered due to remaining 
uncertainties in the MOS response below this energy. 



3. Results 

3.1. Overall spectrum 

The overall MOS1 and MOS2 spectra, within 10' of the 
cluster center, excluding the central R < 2' cooling flow 
region (see below) is shown in Fig. ||. The spectra are cor- 
rected for vignetting and background subtracted as de- 
scribed above. They are compared with the PSPC spec- 
trum of the same region. As the PSPC FOV is much larger, 
the background spectrum can be reliably estimated from 
the outer region (25' — 42'), free of cluster emission. The 
XMM-Newton and PSPC data are jointly fitted in the 
[0.3 — 10] keV range. We let the relative normalization 
between the various instruments be free but assumed a 
common temperature and abundance. When the Nn value 
is frozen to the 21 cm v alue ( A^h = 1-05 x 10 20 cm -2 ; 
Mittaz, Lieu & Lockman 1998 ), the best fit gives kT = 
5.95 ± 0.1 keV and an abundance of 0.27 ± 0.02. The re- 
duced x 2 is 1-27. If the Nh value is allowed to be free, 
we derive Nh — 0.95 ± 0.07 x 10 20 cm~ 2 , in agreement 
with the 21 cm value and the x 2 value, temperature and 
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Fig. 5. Comparison of XMM (open circles) and PSPC (full 
line) surface brightness profiles. The profiles are background 
subtracted and corrected for vignetting (see text for details). 
The PSPC data are renormalized to the total number of XMM 
counts in the 2' to 10' region. Insert: zoom of the region where 
the discrepancy is maximum. Note the excellent agreement in 
shape. 

abundance are unchanged. In the following we thus fix the 
TVh value to the 21 cm value. 

The reduced \ 2 is reasonable, the deviations between 
model and data are less than typically 10%^]. One should 
also note the excellent consistency between the MOS1 
and MOS2 spectra (relative normalization of 1.02) and 
between the PSPC and XMM (relative normalization of 
1.05) The best fit overall temperature in this region is in 
good agreement with the global temperature value from 
ASCA (6.0 ± 0.3 keV, Markevitc h et a l. |l998j ) and SAX 
(6.0 ± 0.4 keV, Irwin & Bregman |2000D . Overall this sug- 
gests that the data processing and the on-axis response 
matrix are reasonable. 

When the best fit model is extrapolated down to 0.1 
keV, the known soft excess in t he EU VE data below 0.2 
keV (Mittaz, Lieu 



Lockman 1998 ) is clearly apparent 
in the PSPC spectrum. This excess cannot be yet studied 
with XMM, due to the previously mentioned calibration 
uncertainties (Section 2.3) in the response matrix in the 
[0.15 - 0.3] keV range. 

3.2. Surface brightness profile 

An independent check of the background subtraction and 
vignetting correction was performed by comparing XMM 
and PSPC surface brightness profiles. The PSPC profile 
in the [0.9 — 2.] keV energy range was extracted and cor- 
rected for vignetting using the standard EXSAS procedure 
for deriving exposure maps. As was done in the spectral 
analysis, the background was estimated from outer regions 
in the FOV with point sources removed. The PSPC pro- 
file is compared in Fig. |^ with the vignetting-corrected 

2 The most obvious misfit is seen around the O edge (E~ 
0.55 keV) in the MOS2 spectrum but not in the MOS1 spec- 
trum. Its origin is unclear. 




Temperature (MOS1) 

Fig. 6. Comparison of the temperature derived from 
EPIC/MOS1 and MOS2 data for the various rings of the ra- 
dial temperature profile. The errors are at the 90% confidence 
level. 



summed MOS1 and MOS2 profile in the [1 - 3] keV band. 
The energy bands were chosen so that the relative emis- 
sivity between XMM and the PSPC were insensitive to 
possible temperature gradients, while at the same time 
retaining a good S/N ratio. The XMM background profile 
was estimated from the Lockman Hole data. The profiles 
are binned so that there is at least a 5 a detection in each 
bin. For direct comparison the PSPC data in the figure 
were renormalized to the total XMM counts over radii 
from 2' to 10'. The agreement in shape is excellent up to 
15' with the maximum difference (at the 25% level only) 
occurring at radii between 8' and 11'. This comparison 
further validates our data analysis method. 

3.3. Radial temperature profile 

Spectra in concentric rings (i.e., annuli) centered on the 
cluster X-ray emission peak were extracted from MOS1 
and MOS2 data. We tried to attain roughly similar preci- 
sion in the temperature estimates from region to region. 
Thus the widths of the various rings were chosen so that at 
least a 5<r detection in the [5 — 10] kev range was reached. 
This was possible for all but the final annulus. In addition 
a minimum width of 30" was set corresponding to the 
90% encircled energy radius of the PSF. Each vignetting- 
corrected and background-subtracted spectrum was fitted 
with an isothermal MEKAL model allowing the normal- 
ization, temperature and abundance to be the only free 
parameters. 

The temperatures determined from the MOS1 and 
MOS2 data are mutually consistent (see Fig. [j|). We thus 
summed the MOS1 and MOS2 data and re-determined the 
temperature in each ring from the combined data. The fits 
are good with reduced x 2 values in the range 0.98 — 1.24. 
The worst % 2 value corresponds to the central rings, where 
we arc affected by the cooling flow. Typical spectra are 



6 



Arnaud et al.: Measuring Cluster Temperature Profiles with XMM/EPIC 




channel energy (keV) 

Fig. 7. Combined EPIC/MOS1&2 spectra of the inner 1' - 
1.5' region (black points), the 2' — 2.5' region just outside the 
cooling flow region (red points) and the outer 8' — 10' ring 
(green points). The full lines are the best fit isothermal model 
for each region. 

shown in Fig. [?]: the inner 1' — 1.5' region, the 2' — 2.5' re- 
gion just outside the cooling flow region and the last ring 
(8' — 10' region). The errors in our temperature estimates 
consist of the quadratic sum of the statistical error (at the 
90% confidence level) and systematic errors estimated by 
varying the background level by ±5%. We also fitted the 
spectra in the restricted 0.3 — 7 keV range, which is less 
sensitive to background and we did not obtain any signifi- 
cant differences in the derived temperatures. We also ten- 
tatively estimate the temperature in an outer region cov- 
ering radii of 10' — 12'. The temperature from this region 
is uncertain and the value should be used with caution. 
The fit is not good (x 2 = 1.8), due to an excess at energies 
above 5 kev and a poor fit at low energies. Both effects 
are likely due to uncertainties in the background (both the 
X-ray background at low energies and CR-induced back- 
ground at high energies) which dominates the emission. 
The error includes systematic errors estimated by varying 
the background level by ±5% and restricting the fit to the 
[0.3 — 5] keV energy range. 

The resulting temperature profile is shown in Fig. || as 
a function of radius both in physical units (upper axis) 
and scaled by the virial radius (bottom axis). The scaled 
radius is defined as r/r 2 oo where r 2 oo = 2.76 Mpc is the 
virial radius for a 6 keV cluster at z = 0.063 (Evrard et 



al. 1996) 



4. Conclusion 

We have determined the temperature profile of A1795 up 
to a typical radius 0.4 times the virial radius. Both vi- 
gnetting effects and background subtraction are reason- 
ably well understood. The temperature profile can be com- 
pared with the profile derived by SA X (Ir win & Bregman 
2000|) and ASCA (Markevitch et al. |l998p . The dramatic 
improvement in accuracy and spatial resolution of the pro- 
files must be noticed. However, our results show that the 
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Fig. 8. Radial temperature profile as a function of angular ra- 
dius (upper axis) and scaled radius (lower axis), derived from 
XMM/EPIC data. The horizontal line corresponds to the mean 
value derived from fitting the overall spectrum of the region 
from 2' to 10' radius. The dotted line is the universal profile 
derived by Markevitch et al. ( 1998| ) from ASCA data, normal- 
ized to the virial radius and mean temperature of A1795. 



limiting factor for determining temperature profiles with 
XMM is the high level of the CR-induced background. 
Improvement of background models would be of great 
help. 

The drop of temperature in the center as measured by 
ROSAT (Briel and Henry 1996), but not seen with SAX, 
is unambiguously confirmed. The drop in temperature is 
also readily apparent in the spectrum of the central region 
shown in Fig.?], which shows the increasing importance of 
the Fc L line complex. Moving in toward the cluster center 
the temperature starts to drop at about 2'; further in the 
temperature profile is, for the first time, well resolved. 
A significant excess in the surface brightness profile (as 
compared to the overall j3 model) is also observed in the 
same region (Briel & Henry 1996). This temperature drop 
indicates that the gas is cooling. A detailed study of this 
inner region, based on both EPIC and RGS data, can be 
found in Tamura et al. (2001). 

Beyond 2' the temperature results are consistent with 
an isothermal profile at the mean temperature (full line in 
Fig.^) in agreement with SAX results. The ASCA temper- 
ature beyond 6' (Markevitch et al. 1998| ) is also in excellent 
agreement with our data, but is significantly higher in the 
1.5-6' radial range (kT = &.2t\% keV). The temperature 
in the central r < 1.5' region is corrected for the presence 
of the cooling flow and cannot be directly compared with 
our data. We also show in Fig.|| the decreasing 'universal' 
profile obtained by Markevitch et al. ( 1998j ) from their 
compilation of ASCA temperature profiles (dotted line). 
It can be safely compared to our data outside the cooling 
flow region. Our results do not support this type of profile 
for A1795. However, we emphasize that no definitive con- 
clusion should be drawn yet, based on the observation of 
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one cluster only and while we arc still in the early phases 
of the XMM mission. In particular we have assumed that 
the (small) energy dependence in the shape of the PSF 
can be neglected. This still needs to be confirmed with 
in-flight calibration data and detailed simulations. 

Acknowledgements. We would like to thank J. Ballet for sup- 
port concerning the SAS software and useful discussions on 
the vignetting correction method. We thank J.-L. Sauvageot 
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